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ABSTRACT 


2 W x 


The shrouded turbine bucket natural frequencies and associated resonant shaft, 
speeds for the M-l Liquid Oxygen Turbopump are presented in this report. Both analy- 
tieal and experimental analysis are discussed and correlated. 

Using the experimental findings, the operating range of 3500 rpm to 4000 rpm 
was investigated for reasonant points. Analysis showed that the ^first-srage bxade 
package could develop a secondary resonant problem as a result of its second mode of 
tangential vibration being excited by the second harmonic of the upstream 43 nozzle 
blade stimulus. No resonant problems are expected for the second-stage blade 
packages , 
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I. 


SUMMARY 


The analytical and experimental determination of the natural frequencies of the 
turbine blades/packages was completed for the first-stage and second-stage rotors of 
the two stage turbine design of the oxidizer turbopump for the M-l Engine, The cor- 
relation reveals quantitive agreement between measured and predicted frequencies for 
the range extending up to 10,000 cps. In addition, mode shape measurements provide 
qualitative identification of the predicted shapes associated with these frequencies. 

Using the natural frequency data for the various blade package combinations, 
the possibility of developing blade operating resonances was studied. Some resonances 
were found in the operating range* They consisted of the higher harmonics of blade 
passing stimulus coinciding with higher mode natural frequencies. These resonances 
are sharp, which means that the stimulus must be very near the natural frequency to 
excite the blades to any appreciable amplitude of vibration. However, the condition 
of narrow resonant response does not entirely eliminate resonances as a potential 
problem. 

II. INTRODUCTION 


The purpose of this analysis was to study the possibility of developing reson- 
ances of the Model II oxidizer turbine blades packages at various shaft speeds. The 
natural frequencies of the blades must be known to predict these resonance points. 
Resonance exits when the fundamental or higher harmonic frequency of the exciting 
force coincides with a natural frequency of the blade. 

Natural frequencies are evaluated analytically by means of a numerical proced- 
ure that performs a detailed analysis of a shrouded group of blades. All of the 
natural frequencies and associated mode shapes (tangential^ axial, and torsional) 
which fall in any specified frequency range can be ascertained. 

For purposes of correlation with the analysis, the spectrum of natural fre- 
quencies for the first-stage and second-stage turbine rotors was determined experi- 
mentally in the range of 50 cps to 10,000 cps. Frequencies and qualitative mode 
shapes were evaluated for the four-blade and five-blade groups in each rotor stage. 
The turbine stator was not tested. 

Ill . VIBRATION ANALYSIS 

A. DESCRIPTION OF THE TURBINE SYSTEM 

1, Actual Turbine Blade Configuration 

The Model II Oxidizer Turbopump Assembly turbine system consists 
of s a first-stage rotor blade, P/N 286528; a stator blade, P/N 286545“ and a second- 
stage rotor blade, P/N 286536. The turbine system layout is shown in Figure No. 1* 
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The Model 21 Oxidizer Turbopump Assembly first-stage rotor blade 
consists of 98 blades banded together in groups of four and five blades. The blades, 
which are made of Inconel ?l8 material are interconnected by a shroud, which Is also 
made of Inconel 718. A detail of the blade and shroud configuration of the first- 
stage rotor blade is shown in Figure No. 2. The rotor blades are fixed to the tab and 
shroud by welding and brazing* Note that the blade section at the hub is reduced at 
the trailing edge. 


The first-stage rotor blade length varies from 3 . 52-in, at the 
leading edge to 3.88-in. at the trailing edge. The shroud is canted at an angle of 
six-degrees, fifteen— minutes. The turbine blade pitch at a mean diameter of 33«00-in. 
is , 997-in. and at the shroud, the pitch is 1.074-in. 

The Model II Oxidizer Turbopump Assembly stator blade, consists 
of 9? blades shrouded together in groups of five and six blades. Because the stator 
blade vibration analysis will not be reported, a detailed description is omitted from 
this report. 


The Model II Oxidizer Turbopump Assembly second-stage rotor blade 
(see Figure No. 3) consists of 94 blades banded together by means of a shroud in groups 
of four and five blades. The blade geometry of the second-stage rotor is identical to 
the first— stage rotor blade. However , the second— stage shroud is slightly shorter 
than the first— stage shroud and it is not canted. The blades are fixed at the snroud 
and hub by means of welding and brazing. Again, a reduction of section properties 
occurs at the hub 0 


The second-stage blades are longer than the first-stage blades. 
They vary in length from 5o5“in<> at the leading edge to 5. 7 -in, at the trailing edge. 
The blade pitch at the mean diameter of 32.8l5“in. is 1.0l6-in. whereas the blade 
pitch at the shroud is 1.196-in. 

2. Simulated Blade Package for the Vibration Test 

A simulated blade package for the first-stage and second-stage 
rotor was manufactured for the vibration tests. The first-stage rotor, P/N 286529, 
(Figure No, 4) and the second-stage rotor, P/N 286332, (Figure No. 5) consist of 
four-blade and five-blade packages banded together by a shroud. The geometry of tta 
blade hub end shroud conforms to the Model XI specifications. 

The blade fixity at the shroud and hub are the same as In the 
design of the Model II blades. For the vibration test, the simulated blade package 
is mounted on a holding fixture (see Figure No. 6) which is then bolted to the MB- 
C60 electro-dynamic exciter. 

B. LUMPED PARAMETER VIBRATION ANALYSIS 
1. Method of Analysis 
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First-Stage Simulated Blade Package, P/N 286529 



Second-Stage Simulated Blade Package, P/N 286532 



Test Set-Up of Simulated Blade Package for Installation on C-60 Exciter 
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The natural frequencies of the banded turbine blades were pre- 
dieted utilizing two Aerojet-General Corporation computer programs * which are for in= 
plane and out-of-plane vibrations of in-plane structures, 

a. Assumptions 

The application of these computer analyses for the fre- 
quency and mode-shape prediction for the banded group of blades requires the following 
assumptions s 

fl) The system consists of a series of buckets, elas- 
tically jointed at the tip by a shroud. 

(2) The principal axes of inertia of all bucket cross- 
sections are oriented in the tangential and axial directions, respectively. 

(3) The center of twist of each bucket cross-section 
coincides with the center of gravity of that cross-section. 

(4) The mass distribution of the buckets can be repre- 
sented sufficiently as lumped masses. 

The lump parameter approach assumes that a mass distribu- 
tion can be idealized by a system of lumped masses at a discrete point. The struc- 
ture to be analyzed is defined as a main system with branches. For example, consider 
the following structures; 
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b. Mathematical Model 

The computer program analyzes a system by relating the 
conditions of state (i.e. Moments, Shear, Axial Forces, Deflections, Rotations, etc.) 
at the first interior point to the boundary conditions at the start; then in turn, re- 
lating the conditions at the second interior point to the first interior point, which 
then relates the second point to the start. This procedure is continued along the 
main system until the joint of the first branch is reached. Then, this relating 
start-to-interior-points technique is used on the branch until the end of Branch 1 
is related to its start. Applying the Branch 1 start boundary conditions at the end- 
to-start relationships, the forces at branch and in terms of its displacements are 
obtained. This allows the conditions of the main system past the branch to be rela- 
ted to the condition before the branch. The procedure is continued until the end of 
the main system is related to the last interior point and hence to the start. By 
utilizing the known boundary conditions at the start and end, the unknown boundary 
conditions can be evaluated. Then, all the interior conditions can be evaluated by 
"re-walking" through the main system and branches. 

To facilitate this "walking through" technique, transfer 
matrices are utilized to get conditions of one end in terms of the other end. To 
illustrate, consider a lumped parameter model broken into bays, a typical bay is shown 
below; 
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Let -CaJ- be the state vector as point n. Then, an elasti- 
city matrix [e] gets the conditions of state at (2) in terms of (J) and the applied 
loading between (2) and (l) . The mass matrix IF] relates (|) to (2) another elas- 
ticity transfer matrix (E] gets conditions at © in terms of Q) , and the angle 
joint transfer matrix [IQ gets (5) in terms of © . Hence, from the above proce- 
dures 


Mjr Mq “ fe] IR ( A j START 

Mi ®(g) “ M x PIjr l A } START 

Mil = W il Mir, iAh start 

= M x M-(g) ■ [ R ] x M il Mi Di^ -[a]- start 

Thus, it follows that for a multiple bay structure, the start to |/\] relationship iss 





Mnl Mh 


“W" START 


After the last bay of the first segment of the main system has been related to the 
start, the first branch can be "walked through" to obtain the branch joint transfer 
Matrix Cb3, thus, crossing over the branch, the state vector becomes, 


■H- v 


X2 


ytli + , ( k w N wj 


H n [A] 


XL 


Continuing on through to the end of the main system, the final state vector becomes ; 


-CaK 


END = 


Vr 3 * 1 

i = 1 


N=,X. 

’/ T 

N = X(i~l)+1 


(&0 N M nl M n [e] nr ) [ b ] i_i 




START 
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Where; [b] q = M unity matrix; [b]. = 1 For ±> 1 ±s as defined in Appendix A. 
j = number of branches 

The state matrix and the transfer matrix derivations are given in Appendix A for in- 
plane and out-of-plane loading. 

c. Use of Symmetry and Anti-Symmetry Model 

When a banded group of blades vibrates in the axial direc- 
tion, the mode shapes tend to take on either symmetric or anti-symmetric appearance. 
Thus, the frequencies for the modes can be computed by considering only one-half of a 
banded group of blades. By using the method of symmetry (i.e. , utilizing appropriate 
boundary conditions at the point of symmetry) , the mode shape of the entire banded 
group can be determined. The above procedure also applies to the blades that behave 
in an anti-symmetric manner. This procedure was used only in the axial vibration 
analysis ; however, it can also be used for tangential vibration analysis. 

A schematic with the boundary conditions for a symmetric 
and anti-symmetric case is shown below; 


V=0 



Symmetry 


r 

i 

i 

i 

i 

rbr 



The anti-symmetry analysis yields the frequencies for the 
odd number modes, while the symmetry method is used to evaluate the even number modes. 


The use of the symmetry and anti-symmetry method of analy- 
sis for axial vibration is limited to a banded group having an even number of blades. 
Thus, for the vibration analysis of the first-stage and second-stage rotor blades, only 
the four-bladed group was determined. However, this method is more accurate than when 
the total system is considered because only half as many computations are required. 


2. Model Representation 

A lump parameter model for the Model II Oxidizer Turbopump Assem- 
bly, turbine bucket blades and shroud attachment are shown in Figure No. 7 for the 


Page 12 




L = bay length 

V - weight 


Figure No 0 7 
Lump Parameter Model 
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tangential and axial vibration analysis,, This model represents a four-blade package 
The model for the axial vibration utilizes the symmetry and anti-symmetry method of 
approach. 


The lump parameters between the first-stage and second-stage 
blades differs not only in the shroud properties but also in the length of the blades 
and the shroud pitch. The length of bays and the blade masses for the computer pro- 
gram are tabulated in Figure No. 7 for both the first-stage and second-stage blades. 

The turbine bucket section properties for both first-stage and 
second-stage are identical. The area and moment of inertia of the bucket sections 
are shown in Table I. The section properties at the bucket fixity to the hub have 
been reduced slightly to simulate the actual section. In essence, this reduction 
lessens the bucket stiffness at the hub. 

The section properties of the shroud for the first-stage and 
second-stage blades are shown in Table II. The shroud configuration used in the com- 
puter program analysis is the actual geometry of the simulated test blade package* 

The section properties include the area, moment of inertias, shear deflection constant 
k, and torsional constant K. 

3. Numerical Results 

a. Mode Shapes and Frequencies 

The mode shapes and frequencies for the Model II Oxidizer 
Turbopump Assembly first-stage and second-stage rotor blades have been determined by 
utilizing the matrix transfer technique and the lump parameter models (see Figure 
No. 7). Further discussion of the tangential, axial, and torsional modes follow. 

(1) Tangential Mode 

The mode shapes and frequencies for the tangential 
vibration analysis were predicted for the first-stage and second-stage rotor blades. 
These are summarized in Table III for the four-blade package and Table IV for the 
five-blade package . 

The results show that between the first and second 
tangential modes, there is a group of frequencies which are called the fundamental 
fixed-support ed modes. These modes, which are the first, second, and third modes, 
are characterized by large vibration amplitudes near the bucket pitch line but very 
little motion at the tips.. 


The numerical results also indicate that the number 
of discrete frequencies in the tangential fixed-support ed mode group is always one less 
than the number of buckets banded together. For example, for a four-blade group there 
are three discrete frequencies. The vibration patterns of the buckets, as observed 
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C = shear deflection coefficient 
k = shear deflection constant 
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TABLE'] III 


CALCULATED TANGENTIAL AND AXIAL NATURAL FREQUENCIES OF FIRST-STAGE TURBINE 


Tangential 
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Axial 
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Axial 
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Four -Bucket Group 
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TABLE IV 


CALCULATED TANGENTIAL AND AXIAL NATURAL FREQUENCIES OF SECOND-STAGE TURBINE 


Tangential 

Mode 

itxial 

Mode 

Four-Bucket Group 

Five-Bucket Group 

Tangential 

Mode Shape 

Axial 

Mode Shape 
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Mode Shape 
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Second Cantilever 


3820 
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First 

Axial 

Group 
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^157 

4257 

1 












Second Cantilever 


















from the mode shapes (see Tables III and IV) take on "odd 11 or "even" symmetry 
appearance depending upon whether corresponding buckets in each half of the group 
are in-phase or out-of-phase. 

(2) Axial Mode 

The mode shapes and frequencies for the axial 
vibration analysis were predicted by using the lump parameter model shown in Figure 
No. 7. This model utilizes boundary conditions consistent with the symmetric and anti- 
symmetric modes. The numerical results of the frequencies for the axial vibration of 
the first-stage and second-stage rotor blades as well as their respective mode shapes 
are shown in Table IV for the four-blade package. 

The first axial group modes are characterized by 
the progressive pattern of nodes in the shroud band. There are no nodes along the 
bucket length. In the lowest or first mode, there are no nodes along the shroud band. 
All of the buckets in the first axial group vibrate in the axial direction and all of 
their motions are in phase. The second axial mode has one node along the band, the 
third has two nodes, and so forth. The highest frequency mode always has nodes in the 
shroud between each bucket in the banded group and thus, the number of nodes along the 
band is equal to one less than the number of buckets banded together. Therefore, the 
total number of modes in this axial group is always equal to the number of buckets 
banded together. 


(3) Torsional Mode 

For shrouded blades pure torsional modes do not 
exist. Moreover, the coupled torsional-bending natural frequencies are as predicted 
for the axial modes of vibration. 

C. VIBRATION TESTS OF FIRST AND SECOND STAGE TURBINE ROTOR BLADES 
1, Test Set-Up and Procedure 

The vibration tests utilizes the simulated blade package (see 
Section III, A) to determine the various resonant frequencies and corresponding mode 
shapes. The first-stage rotor, P/M 286529* and second -stage rotor, P/N 286532, each 
in turn, was mounted in a holding fixture and bolted to the MB-C60 electro-dynamics 
-;>:citer (see Figure No. 6). The simulated blade packages were excited in both the 
axial and tangential directions from 100 cps to 9000 cps. 

The resonant frequencies were determined by the phase change 
between the accelerometer output and the excitation input as observed on an oscillo- 
scope. An Endevco accelerometer was used to pick up the output signal on the blades. 

The nodal patterns of the banded blade were determined by using 
the salt pattern technique and the damping-undamping technique. 
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Test Results 


The mode shapes and resonant frequencies from the vibration tests 
of the first-stage and second-stage turbine blades are summarized in Tables V and VI. 
For the tangential modes, the first cantilever, first fixed-support ed, and second 
cantilever modes were obtained for the first-stage and second-stage turbine blades. 

For the axial modes, only the first axial group was measured from the vibration test. 

The frequency values in the first fixed-supported group for both 
turbine blades were quite close together, whereas in the first axial group, the fre- 
quency varied widely, yet they are lower than the first fixed-supported group. 

The test results also indicated that the frequencies of the four- 
blade and five-blade packages in the higher modes are similar with the exception that 
the five-blade package has one more mode than the four-blade package. Thus, use of 
the four-blade package data is sufficient for evaluation of resonance. 

3. Correlation of Test Frequencies with Analytical Predictions 

A comparison between calculated and observed frequencies of 
banded bucket groups is presented in Table VII for the four-blade and five-blade 
banded groups of the first-stage and second-stage rotor blades. 

In the tangential mode, comparison of frequencies is given for 
the four-blade and five-blade groups, whereas in the axial direction only frequencies 
from the four-blade banded group were measured. The results, which are very encourag- 
ing, indicate that fairly accurate predictions of resonance conditions can be accom- 
plished through the use of these digital computer programs. Therefore, these computer 
programs can aid in the future design efforts of banded buckets. 

The discrepancies between the test values and the computer pre- 
dictions in the higher modes could largely be the result of the blade orientation? 
that is, the shroud is actually of such a nature as to couple, although weakly, the 
tangential and axial modes of vibration. Also, the treatment of the shroud stiffness 
and "effective" length becomes more significant in the prediction of blade frequencies 
in the higher harmonics. In the analytical prediction analyses, the effective length 
influence was accounted for by distributing the shroud stiffness as shown in Figure 
No. 8. 

D. RESONANT FREQUENCY DETERMINATION 
1. Excitation Frequencies 

a. Description of Primary Stimulus for Each Stage 

The natural frequencies of the turbine blades are excited 
by the action of periodic stimulus which, in this case, is a gas flowing through a 
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TABLE V 


OBSE RVED TANGENTIAL AND AXIAL NATURAL FREQU ENCIES OF FIRST-STAGE TURBINE 
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TABLE VI 


OBSERVED TANGENTIAL AND AXIAL NATURAL FREQUENCIES OF SECOND-STAGE TURBINE 
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TABLE VII 


COMPARISON OF OBSERVED AND CALCULATED TANGENTIAL AND AXIAL NATURAL FREQUENCIES 
OF FIRST-STAGE AND SECOND-STAGE TURBINE BLADE 
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group of nozzles. The stimulus has a period of time equal to the time required for 
the bucket to travel a distance of one nozzle pitch. Because the stimulus is non- 
sinusoidal, it can also excite the buckets at frequencies that are multiples of the 
frequency of passing nozzles. 

The frequency of the stimulus may be calculated by 

n NG 

Frequency s — ’ , cps 

where s n = the harmonic 

N = number of nozzle per 360 degrees 
SI = shaft speed, rpm 

Now, by making a plot of frequency, Q , (cps) , versus shaft speed, (rpm), a series of 
lines representing the frequencies of the stimulus for each harmonic can be drawn, 
and the slope of these lines is defined ass 


Freq of Stimulus n N 160) - n N 

Slope - Shaft Speed Si 60 

The primary stimulus for the first-stage turbine blades is the result of gas flowing 
from the 43 turbine nozzle blades. The second-stage turbine blades can be excited by 
the actions of both the 97 stator blades, which is the primary stimulus, and the 43 
nozzle blades (see Figure No 0 l) 0 

Plots of the stimulus for the first-stage and second-stage 
turbine blades are shown in Figures No. 9 and No. 10. Resonant points of the turbine 
blades (i.e. , the point in which the natural frequency coincides with the frequency of 
the stimulus) are determined by the intersection of the horizontal lines representing 
the natural frequency of the blades with the stimulus lines. These maps represent 
"critical speed" plots for the turbine blades. The horizontal lines would become 
curved lines if centrifugal stiffening of the blades is significant! however, this 
effect was determined and found to be negligible. 

2. Tangential and Axial Mode Resonance Points 

Critical speed plots for the first-stage and second-stage turbine 
rotor blades are given in Figures No. 9 and No. 10. Because the operating range of the 
shaft speed is 3500 rpm to 4000 rpm, only the resonance points within the operating 
band are of interest. It should be noted that for each case, the stimulus frequency 
curves are for the first three harmonics (n = 1,2,3). The amplitudes for the higher 
harmonics should be very small and can be considered negligible. The turbine blade 
natural frequencies for the tangential and axial modes are the measured values (see 
Tables VI and VII) of both the four-blade and five-blade groups. The measured 
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frequencies have been corrected to the operating condition by the ratio of the 
elastic modulus at operating temperature. Thus, at a temperature of 1100 F, the 
frequencies are corrected by (.945) x room temperature frequency. 

For the first-stage turbine blades , the critical speed plot . 
(see Figure No. 9) indicates that the second harmonic of the exciting force provided 
by the 43 nozzle blades could theoretically excite the second tangential mode of 
vibration. The natural frequency for this case is 5270 cps, which corresponds to 
the first mode in the fixed-support ed group. 

The critical speed plot of the second-stage turbine blade 
(Figure No. 10) indicates several resonance points within the operating range, par- 
ticularly the first harmonic of the secondary stimulus from 43 nozzle blades. The 
primary stimulus is provided by the 97 reversing blades. The first harmonic of t e 
exciting force provided by the nozzle blades could theoretically excite the first 
fixed-supported group of the tangential modes and the third mode of the axial vibra- 
tions in the five-blade package as well as the similar modes of vibration in the 
four-blade package. At a higher frequency, the third harmonic of the exciting 
force from the 43 nozzles could excite the fifth mode (7900 cps) of the axial modes 
of vibration in the five-blade package. These higher modes of possible resonance 
are not considered critical because they are more difficult to excite and the ex- 
citing stimulus magnitude decreases rapidly with the higher harmonics of the forc- 
ing frequency. 

3. Operating Speed Range Potential Vibration Problems 

The operating range of the Model II Oxidizer Turbopump Assem- 
bly turb in e blade is between 3500 rpm and 4000 rpm, Within this operating range, 
several resonant points occur, particularly with the second-stage turbine blade. 
However, the excitation force of the second-stage blade has been provided mostly 
from the 43 nozzle blades which, in this case, is not the primary stimulus. There- 
fore, the excitation force is not very powerful and hence, a resonance problem is 
not expected,. 

For the first-stage turbine, the resonant points within the 
operating range consist of the second and third harmonic of nozzle passing stimu- 
lus exciting the second and sixth modes of tangential vibration. As no fundamen- 
tal modes have resonant points in the operating range, the turbine blades can 
operate within the range without any significant resonant problems. If a problem 
does occur, it probably will be associated with the second mode resonant point 
excited by the second harmonic of nozzle passing stimulus, 

IV. CONCLUSIONS 

The analytical and experimental determination of the natural frequencies of 
the first-stage and second-stage turbine blades of the Model II Oxidizer Turbopump 
Assembly has been fully evaluated. The results from the vibration tests of a 



simulated blade package, as compared with the analytical solution, demonstrate that 
the natural frequencies from a shrouded group of blades can be predicted quite well 
by the technique presented herein. 

Discrepancies between the predicted and test data could possibly be the re- 
sult of the principal bending axis of the blade package not being in the same plane 
as that of the shroud. For simplicity, in the computer solution, the principal bend- 
ing direction for both the blade and shroud is taken to lie in the same plane. 

Good correlations were obtained in the lower frequencies for both the first- 
stage and second-stage turbine blades. For the higher modes, natural frequencies 
of axial vibration predictions tend to be somewhat higher than the test results. 

The resonance points of the turbine blades using test results were investigated 
in the operating range of 3500 rpm to 4000 rpra. For the first-stage blade package, . a 
resonance problem might occur because of the second mode of tangential vibration being 
excited by the second harmonic of the 4-3 nozzle blade stimulus. Moreover, this is no 
considered to be a primary resonant problem as it is not definite that sufficient ex- 
citation will occur. No resonant problems are expected for the second-stage blade 
packages . 
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APPENDICES 




APPENDIX A 


NOMENCLATURE 





A 

C 

G 

S 

I 

J 

K 

k 

L 

M 

Q 

T 

u 

V 

y 

w 

o4 

$ 

© 

CO 

[B] 

m 

m 

&] 

[A] 

AE 

BG 

El 


2 

= Area , in 

2 

= Shear deflection coefficient, l/in 
= Modulus of Rigidity, lb/in^ 

2 

= Acceleration of gravity, 386 in/sec 

4 

= Moment of Inertia, in 

2 

- Rotary inertia constant, 91 'bs-sec -in. 

= Torsional constant 

= Shear deflection coefficient 
= Length of bay element, in, 

= Bending moment, in-lbs 

- Torque , in-lbs 

= Axial force, lbs 

= Displacement perpendicular to axis of bay element, in. 
= Shear force, lbs 

= Displacement parallel to axis of bay element, in. 

- Weight of bay element , lbs 

ss Pitch angle between bay elements, degrees 
= Rotation, radians 

= Torsional Rotation of axis of bay element, radians 

- Frequency , cps 

= Branch Transfer Matrix 

- Elasticity Transfer Matrix 
s Mass Transfer Matrix 

as Joint Transfer Matrix 
ss State Vector 

2 2 

ss Beam Axial Stiffness Parameter (in )(lb/in ) 

b .2 

& Torsional Stiffness Parameter, (in )(lb/in ) 

2 if 

ss Beam Stiffness Parameter, (lb/in ) (in ) 


Subscripts 


A 

L 

N 

R 

T 


-- axial direction 
= left side of bay element 
a element number 
= right side of bay element 
=s tangential direction 
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APPENDIX B 


VIBRATION ANALYSIS USING 
MATRIX TRANSFER METHODS 






I. VIBRATION A NALY SIS OF IN-PLANE LOADING OF IN-PLANE STRUCTURES 

WITH BRANCH 

The general procedure of the matrix transfer method of analysis is 
outlined on pages B-5 through B-8 of this appendix,. The detailed derivations 
and definitions of the elasticity, mass, and branch joint transfer matrices 
are also set forth* 

The equations of equilibrium and compatibility for the elasticity 
element and lumped masses as well as the joint transfer equations for in- 
plane loading are as follows. 

A. ELASTICITY TRANSFER MATRIX 


V 

NR 
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In matrix form, the elasticity transfer matrix "becomes; 



B. MASS TRANSFER MATRIX 

v 

~K~ 



U NL = U NR’ $RL = $NR’ Y NL = Y NR 
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In matrix form: 



0r { A }hL = [ F ] N Wire 


C. JOINT TRANSFER MATRIX 
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In matrix forms 



° r { A }nL t R ] H |a]-hr 


D. BRANCH JOINT TRANSFER MATRIX 



Starting at Point A "walking" to B with the aid of the 
transfer matrices we get 

B = ^NL Mnr) A Equation (1) 

for R Bays in the branch. 
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Expanding Equation (l) 





K Pf {cQ. 

•N »N 

< S3 . < p- 

< S < w> 

•p 

CD 

PI 


n 



Eh 




XA 

I 

PQ 


0 ) 


SP 


Px 



Expanding Equation (U) gives 


p 

- Hi 

K 

u 

+ Ha 

K 

g 

+ d 1T 

A 

w 

= £ 21 

K 

u 

+ Fj 

^22 

K 

g 

+ dg^ 


Equation (5a) 
Equation (5b) 


-1 


By multiplying Equation (5b) by £ we get 

■ 21 


£/" *w = K +£ 1 £.K +6 .d 

21 U 21 22 S 21 27 


or 


>. -1 . -1 . >• -1 
K = C w - £ £ • K - c • a. 

u 21 21 22 & 21 


‘27 


Equation (6) 


Substituting Equation (7) Into Equation (6a) and writing in the partitioned 


matrix form, we haves 


l t ' 1 1 £ - £ l 

11 21 , 12 11 21 W 22 



-1 


17 


-e. e 


11 


Expanding Equation (7) gives 




✓s. 
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A 

M 


- 


V 1 I B 


✓A 

11 

*12 

H.3 | 

A 

k lU 

A 

k 15 

A 

k l6 

1 k 17 

Hi 

Ha 

k 23 I 

| \h 

A 

k 25 

k 26 

1 k l8 

k 31 

A 

k 32 

A 

V 

33 


k 35 

k 3 6 

1 A 

1 k l9 

. 




— 

— 

. 

4 - - 

0 

0 

0 

I 0 

0 

0 

i i 


-1 

21 


27 



_w_|| Equation 

(?) 


B 

A 

a B 
<*T| 

* /v 


Equation (7a) 
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Hence, we now have the forces at Point B in terms of the displacements at B 
and the known boundary conditions at Point Ao However, these values are with 
respect to the sign convention of the branch; therefore, a change of coordinates 
transfer matric is needed. For simplicity, rewriting Equation (7) and 
Equation (7a) as 


P A 


K 11 

K 12 

k 13 


w 1 

__B_ 



> 




j 

K 

g 

i Equation (8) 

i 


0 

0 

1 


___ 


\ 4 






J 



and multiplying by the appropriate coordinate transfer matrix J yields 


-1 -1 , -1 
G 1 ^ 11^21 G 1 ^ G 1^12 " G 1 ^ 11^21 ^22 G 1 d 12 ” G 1 ^ 11 ^ 21 : d 27 


/ 1 

r 



B 


v 

>=■■■ 

1 



L. 


0 


— t 


0 


— 



(* s \ 


W B 

‘ j 



1 

_ 

J 


\ Equation ( 




cos T 

sin Y 

0 

where : 

• G 1 " G 2 “ 

-sin Y 

cos Y 

0 



0 

0 

1 
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Expanding Equation (9) we have imposed on the main system from the branch 


V 


k ll 

k 12 

k l3 

k i4 

k 15 

he 

H 

-<3 

J 


'V 

Ub 


V B 

< 

) = 

k 21 

k 22 

K 23 

k 24 

k 25 

k 26 

k 27 

< 

$B 

«? 

) Equation (10) 

m b 


k 31 

k 32 

k 33 

k 34 

k 35 

k 36 

k 

37 




1 

v ) 


0 

0 

0 

0 

0 

0 

1 


Q f 

> 



We are now ready to write down the branch joint transfer matrix- Consider the joint 
point of the main system. 



The equilibrium equations are: 

V iL * V iH * V *» •- *tt + V M i X “ M 1R + “b Ration (U) 
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branch joint transfer matrix is? 
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II. 


VIBRATION ANALYSIS OF OUT-OF-PHASE LOADING OF IN-PLANE 
STRUCTURES WITH BRANCHES 

The equations of equilibrium and compatibility for the elasticity 
element and lumped masses and the joint transfer equations for eut-of -plane 
loading are as follows. 

A. ELASTICITY TRANSFER MATRIX 
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In Matrix Form: 
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o 
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D. BRANCH JOINT TRANSFER MATRIX 

The branch joint transfer matrix equations for out- 
of -plane leading are derived in the same manner as for the in-plane 
loading; therefore, refer t© section I, D, page B-Uo 
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